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ANALYSIS OF A RADIAL-FLOW HALL CURRENT 
MAGNETOHYDRODYNAMIC GENERATOR 


by Lester D. Nichols 
Lewis Research Center 

SUMMARY 

A magnetohydrodynamic generator of two-disk geometry operating in the Hall mode 
is considered. The generator working fluid flows radially between the two disks and a 
constant magnetic field is applied parallel to the axis of the disks. The conditions where 
the Mach number is small compared to one and where viscous effects can be neglected 
are analyzed. The power output, power-output density, and efficiency of this generator 
are determined and compared to the linear Hall generator. The calculations indicate that 
the linear generator has higher efficiency and power output for a given working fluid and 
entrance velocity. 


INTRODUCTION 

A magneto hydrodynamic (MHD) generator operating in the Hall mode (i. e. , the Hall 
current is the load current) is considered. Generators of this type are particularly well 
suited for operation with large Hall parameters (ref. 1). A Hall generator with the Far- 
aday electric field equal to zero can be shown to give a large amount of nonequilibrium 
ionization (ref. 2, eq. (23)). In linear generators the Faraday field is made to vanish by 
segmenting the electrodes and then shorting opposite pairs of segments. The short cir- 
cuit for the Faraday current may lead to difficulties (ref. 3) so that a geometry having 
the Faraday field always equal to zero is of interest. Such a configuration is shown in 
figure 1 where the working fluid flows radially between the two disks with a constant 
magnetic field applied parallel to the disk axis. The design is functionally similar to 
that reported in reference 3. 

This report analytically compares the performance of a disk and a linear Hall gen- 
erator in terms of power density and efficiency. Emphasis is placed upon parametric 
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Figure 1. - Sketch of radial-flow generator. 


performance evaluation that illustrates the effect of the 
the geometrical arrangement. Several simplifying as- 
sumptions are used; namely, that viscous effects are 
negligible, magnetic Reynolds number and Mach number 
are small, and electrical properties are independent of 
magnetic field and current. Although these assumptions 
may not always be valid (especially if electron heating 
is imparted), they are made for both generators so that 
the analysis is selfconsistent and should permit a rea- 
sonable comparison of the configuration effects. 


EQUATIONS DESCRIBING GENERATOR 

The momentum equation for an inviscid fluid carrying a current and flowing in a 
magnetic field is (ref. 4) 


P — + grad p = j x B (1) 

Dt 

(a list of symbols is given in appendix A). If the magnetic Reynolds number is small 
compared to 1 then the induced magnetic fields can be neglected in comparison to the 
applied magnetic field strength B (ref. 4). For this study B will be only in the 
z direction. (See fig. 1 for coordinate sketch.) The velocity and current density vectors 
can be written: 

u = v? + w cp 



so that equation (1) in component form becomes: 

( 2 ) 
( 3 ) 

For a flow with the Mach number small compared to one the continuity equation for 
axially symmetric flow requires: 
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P [v — - — 

dr r / dr 


+ * = j B 

J cp 


pv- — (rw) = -j B 

r Hr 1 
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A (rv) = 0 
dr 


( 4 ) 


For a plasma with nonzero Hall parameters for ions and electrons and with the cur- 
rent flowing in a plane perpendicular to the magnetic field, Ohm's law becomes (ref. 5, 
eq. (49)): 


(1 + Vi) 


cr 0 (E + u X B) + /9 e 


B xj 
B 


This equation can be solved for j (ref. 5, eq. 50): 


j = O’ 


(1 + /3 /3 )(E + u X B) + B B X (E + u X B) 
“ a c n 


where 


d + ^i > 2+(3 e 


For the coordinate system used herein, the equation for j can be written in component 
form as shown: 


Jr = a 

[(1 + Vi )(E r + wB) + /3 e vB ] 

(5) 

V = CT 

(1 + ^^(-vB) + /3 e (E r + wB)] 

(6) 

The requirement that div j = 0 is 



— (rjj = 0 

dr r 



which yields 

rj r = (7) 

This value for rj r may be substituted into equation (5) which, when solved for E r 
yields 
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( 8 ) 


r l J l 

1 - - al 3gVB 

E = — wB 

r o(l + W 

Using this value for E r the azimuthal current j may be obtained from equa- 
tion (6) 
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Pjjr - V 6 
l + ^i 


Equations (4) and (3) can be solved for v and w 


v = 


V 1 



r 



p v i 



( 9 ) 


( 10 ) 


( 11 ) 


Equations (7), (8), (9), (10), and (11) express all components of current, electronic 
field, and velocity as functions of radius in terms of initial conditions at r = r^. 


POWER-OUTPUT DENSITY 


The equations in the preceding section describe a generator (positive power output), 
only for a certain range of initial current density, jj. The smallest of interest is 
the open-circuit value, namely zero. The largest value is the short-circuit current, 
that is, the current which flows when the load voltage is zero. The load voltage can be 
expressed as 


V = 



(-E r )dr = 


"l 1 " FJ Vl B '' 


1 + ^ 
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+ + W 5f ] 




( 12 ) 


where 

4 


■ llllll I iMHIHIIiail I 


I I I 



2 



CT o B2r i 

O = ~ 

2^1 + W 2 + <£] 

It is evident from equation (12) that the short-circuit current is 

oViB/3 

Q ) — — ® 

1 sc 1 + 6f(l + (3^) 


It is convenient to express the load voltage as a fraction of the open-circuit voltage. Let 
this fraction be called the load parameter and be denoted by K, so that from equa- 
tion (12) 


K s 



[l + 5f(l + 0 e ^i>] 

°Vi B ^ e 


( 13 ) 


Thus the range on can be specified in terms of the load parameter as 

oVjB^l - K) 


h = 


0 < K < 1 


1 + 6f(l + /3 e /3.) 

The power output can be represented as the volume integral of [j • (-E)] 

P = Xol [1 • (- 5 >] dT - fv[ 2 ( 2 't>i r )(-E r )r dr 


(14) 


With equations (7), (12), and (14), this power becomes 
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P = 27rhr^jjV = 
Since the volume of the generator is 


77hr 


2„2 0 2 t 


\ Inf— J ctv|B p^K(l - K) 


(1 + £ e ^j) [l + 6f(l + ^ e /3.)] 

77h(r| - r 2 y the power density becomes 


n = 



crv 2 B 2 /3 2 K(1 - K) 

(l + ^p^ + sfd + ^.jjd + f) 


(15) 


The maximum value for n occurs for K= 1/2 for the case of constant conductivity. 

The effect of this constant conductivity restriction can be illustrated by comparison 
with results in reference 5. In figure 5 of that reference, the power density for a gen- 
erator with a plasma whose properties are calculated with the theory of electron heating 
is shown to be a maximum for load parameter of about 0. 42. This difference must be 
noted when the generator is operating under conditions such that electron heating is 
likely to be important. 


EFFICIENCY 

The efficiency of this generator may be defined as the ratio of the output power to 
the flow work per unit time done by the fluid. Since the only force considered herein is 
the electromagnetic -body force, the flow work per unit time is j^-u • (T X B)] so that 
the efficiency rj becomes 


f T • (-E)dr 

-'vol 

f [-U • (T x B)]di 

*v Ol 


77 0 K(1 - K) 


(e - A)K + AK + 1 


(16) 


where 


^0 


# + 6 fd + /y*.)] 


(1 + P&) 2 [l + 6f(l + 0/.)] 2 + [/3 e 6f(l + /3 e /3.) 
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1 + 3(1 + p B .)6f 
l + U + jS^M 0 

X + 7] 0 1 + 2(1 + P P.)61 

€ = — = Vo 

2 1 + (1 + P e Pj)6f U 

Differentiation of equation (16) shows that the efficiency is a maximum when 


Vl + e + 1 


( 17 ) 


and has a value given by 


^max 


% 

X + 2( Vl + e + 1) 


(18) 


SPECIFICATION OF MAGNETIC FIELD 

Although the parameters appearing in the expressions for the generator are arbi- 
trary, there are some optimum values that can be specified. Because of ion slip there 
is a magnetic field which maximizes the efficiency and power density (ref. 6, figs. 5, 

6, and 7). To examine this, introduce the following dimensionless variables: 

ff n r i 

6 = (19a) 

2 Vi pv l 

JU 

ju = — (19b) 

Mi 

z = ^ i = (M e B)(M i B) (19c) 

where 6 C is a modified interaction parameter; /i is the ratio of electron to ion mobil- 
ity; and z is an ion slip term. The parameters 8 C and p. do not depend upon B 
since the plasma properties M e > Mp and are assumed constant. Then by substitu- 
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tion in equation (15) the power density may be written as 


2K(1 - K)pv?p6 z 2 

n = 1 c . - ( 20 ) 

rj(l + z)(l + f) [(1 + z) 2 + (p + 6 c f)z + 6 c fz 2 ] 

The value of z for which II as well as the total power output is a maximum is deter- 
mined from 0n/3z = 0 


6 c fz(l + z)(l - z) + pz = (z + l) 2 (z - 2) (21) 

The power output of the generator (but not the power density) also has a maximum 
with respect to variation in the radius ratio. This maximum occurs when 


6 c [(y - 1) - 2f]z(z + 1) - pz = (1 + z) 2 


( 22 ) 


where 


y = 



Equations (21) and (22) determine conditions for which the power output is maximized. 
Thus, for specified 6 C and p (which are parameters determined by the fluid proper- 
ties, generator inlet size, and generator operating conditions) the optimum magnetic 
field and outside radius may be obtained from 


6 c = 


(z 2 - 1) 


z [y - 1 - (z + l)f] 


ix = 


(z - l)(y - 1 - 2f) _ j 
y - 1 - (z + l)f 


(z + D' 

z 


(23) 


(24) 


These functions are shown in figure 2. The resulting maximum power output (divided by 
the input kinetic -energy rate) is given by 
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Figure 2. - Radius ratio as function of modified interaction parameter for specified mobility ratio and magnetic 
field parameter. 


P 5jtzz In y I 

27 rhr jPv j 4(1 + z) [l + (1 + z)0f] 

and is shown in figure 3. The maximum efficiency, calculated from equation (18), is 
shown in figure 4. The ranges shown for the dimensionless variables are those that 
might occur in actual generators. 

COMPARISON WITH LINEAR HALL GENERATOR 

The power-output density for a linear Hall generator is given as (ref. 7, 
eq. (14.43)) 


Mobility ratio,— 
M ~ 
100 


300 

1000 



10"1 10 ° 10 1 10 2 
Modified interaction parameter, 6 C = ao r l^WP v l 


Figure 3. - Power output as function of modified interaction parameter for specified values of mobility ratio. 




<**<1 - K f )vf^B 2 


1 + ^ 


(26) 


It can be shown from equations (20) and (26) that, since 5 c f > 0 and the 

power density ratio Il/IM < 1. This ratio is plotted in figure 5 using the results of 
equations (23) and (24). It may be seen that the ratio decreases with increasing values 
of the slip parameter. 

The efficiency of the linear generator is (ref. 7, eq. (14. 50)) 
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Power density ratio, TT/TT j 



Modified interaction parameter, 6 C = ^± 12 ^^ 

Figure 4. - Efficiency as function of modified interaction parameter for various mobility ratios. 



Slip parameter, (3 e , 


Figure 5, - Power density ratio as function of slip parameter 
for various mobility ratios at optimum magnetic field and 
radius ratio. 



This is equal to the disk- generator effi- 
ciency when f = 0, and has a maximum 
when the load parameter is 



The efficiency at this value of be- 
comes 
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The ratio of efficiencies can be expressed as 




Slip parameter, p e ,Pj 


Figure 6. - Efficiency ratio as function of 
slip parameter for various mobility 
ratios at optimum magnetic field and 
radius ratio. 


and, since from equation (17) 

g 

77 n < e < X and e < 

d + (Vi> 2 

then 

JL< 1 

This ratio is plotted in figure 6, using the values 
from equations (23) and (24). The efficiency ratio 
is seen to be an increasing function of the slip 
parameter . 


CONCLUDING REMARKS 


The results of the calculations indicate that the disk generator is quite inferior in 
performance to the linear generator. For small values of the slip parameter the effi- 
ciency is much less and for large values of the slip parameter the power density is much 
less. However, it must be remembered that this comparison was made on the basis of 
maximum power output of the disk generator. Nevertheless, it appears that the disk- 
generator performance is so inferior that the greater benefits of the electron heating (as 
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indicated in reference 3, (figs. 6 and 10) are inadequate to make its performance com- 
petitive with that of the linear generator. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 14, 1965. 
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APPENDIX - SYMBOLS 


magnetic field strength electron and ion Hall parameters 

electric field strength 6, 6 C interaction parameters, eq. (19a) 

f function of radius ratio e,A, ?7 q efficiency parameters, eq. (17) 


h 

spacing between disks 

V 

efficiency 

r 

current density 

M 

mobility ratios, eq. (19b) 

K 

load parameter 

M e ,Mi 

electron and ion mobilities 

p 

power output 

n 

power densities 

p 

pressure 

p 

density 

r 

radial coordinate 

a ’°0 

electrical conductivities 


unit vectors in coordinate 

dr 

volume element 


directions 

Subscripts: 

u 

velocity 

£ 

linear generator 

V 

load voltage 

oc 

open circuit 

V 

radial velocity 

sc 

short circuit 

w 

azimuthal velocity 

r 

radial component 

y 

2 

square of radius ratio (r/rj) 

1 

inlet condition 

z 

magnetic -field parameter, 
eq. (19c) 

<P 

azimuthal component 
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